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I.  SUMMARY 


The  contribution  of  the  spall  source  to  the  short  period  regional  explosion 
seismogram  is  studied  using  new  parametric  models  for  the  spallation  seismic 
source  derived  from  two-dimensional  (2-D)  axisymmetric  nonlinear  finite 
difference  simulations  for  contained  underground  nuclear  explosions.  It  is  con¬ 
cluded  that  for  high-velocity  crustal  structures  spall  may  be  a  significant  source 
of  short  period  Lg  and  may  significantly  affect  the  Pg  excitation. 

Two  seismic  source  representations  for  spall  are  compared.  Regional 
phase  excitation  from  a  vertical  point  force  and  a  horizontal  tension  crack  are 
examined.  A  new  general  derivation  for  the  equivalence  between  the  two 
representations  at  low  frequencies  (horizontal  wavelength  »  spall  depth)  is 
given.  Numerical  calculations  show  that  the  vertical  point  force  at  the  free  sur¬ 
face  above  the  spall  is  a  good  approximation  in  the  short  period  bandwidth  but 
degrades  at  high  frequencies  (horizontal  wavelength  <  half  the  spall  depth). 

A  smoothing  operator  is  presented  to  account  for  the  distributed  finite 
nature  of  the  spall.  This  operator  has  the  desirable  characteristics  that  it 
represents  a  distributed  spall  disk  with  zero  motion  at  the  tip,  has  an  upward 
concave  spall  surface.  It  falls-off  proportional  to  ar5/2  at  high  frequencies. 
Since  the  spall  source  is  proportional  to  co2  for  low  frequencies,  the  far-field  sig¬ 
nal  should  be  signficant  only  for  a  narrow  band  of  frequencies  compared  to  the 
broadband  explosion  signal.  This  narrow  band  can  be  expected  to  coincide 
with  the  traditional  short  period  seismic  band  (0.5-2  Hz)  for  explosions  in  the  1 
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to  150  Kt  yield  ra.  ge  at  normal  depths  of  burial. 

The  parameterized  model  of  Barker  and  Day  (1990)  is  then  used  to  exam¬ 
ine  the  significance  of  spall  signature  in  regional  waveforms  compared  to  the 
explosion  contribution  to  the  seismogram.  Green’s  functions  are  computed  for 
a  high-velocity  near-surface  layer  crustal  model.  It  is  found  that  for  reasonable 
amounts  of  spall,  the  short  period  Lg  signal  is  completely  dominated  by  the 
spall  source.  The  Barker  and  Day  spall  source  for  a  125  Kt  explosion  buried  at 
680  meter  depth  corresponds  to  a  6.2  mb(Lg )  in  the  Stevens  (1986)  crustal 
model  for  Eastern  Kazakhstan.  The  Pg  signal  may  be  dominated  by  the  spall 
source  for  a  narrower  bandwidth.  Simulations  suggest  that  if  spall  is  a 
significant  source  of  Lg  excitation  then  mb(Lg)  should  be  sensitive  to  scaled 
depth  of  burial. 
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II.  INTRODUCTION 


It  was  shown  by  Day,  et  al.  (1983)  that  a  proper  seismic  source  represen¬ 
tation  for  spall  must  conserve  momentum  and  therefore  the  spall  signal  has  no 
"zero"  frequency  source  strength.  Day,  et  al.  then  showed  that  for  low  fre¬ 
quency  surface  waves  an  equivalent  force  model  for  spall  could  be  represented 
by  a  vertical  point  force.  They  demonstrated  that  in  the  low-frequency  limit,  the 
vertical  opening  of  a  horizontal  tension  crack  is  equivalent  to  a  vertical  force  on 
the  surface  with  a  suitable  time  history.  In  recent  years  this  spall  model  has 
been  increasingly  used  to  represent  the  seismic  spall  signal  at  higher  frequen¬ 
cies  than  it  was  initially  applied.  Stump  (1984)  has  shown  that  for  the  purposes 
of  source  inversion,  suitable  combinations  of  vertical  point  forces  and  explosion 
sources  can  not  be  distinguished  from  an  explosion  plus  a  tension  crack.  In 
subsequent  work,  Stump  (1985)  has  used  the  Day,  et  al.  model  combined  with 
an  ad-hoc  smoothing  operator  to  model  small  contained  chemical  explosions. 
Recently,  Patton  (1988)  has  used  the  Day,  et  al.  model  as  modified  by  Stump 
to  model  regional  Rayleigh  and  Lg  waves  by  moment  tensor  inversion  for  the 
NTS  explosion  HARZER.  Patton  concluded  that  the  spall  signal  was  a 
significant  contributor  to  the  Lg  signal  below  1  Hz.  Johnson  (1988)  using  near¬ 
source  data  performed  moment  tensor  inversions  for  HARZER  and  another 
NTS  explosion,  CHANCELLOR,  and  found  similar  results  to  Patton.  In  both 
cases,  it  was  found  that  a  source  secondary  to  the  explosion  could  be 
represented  as  either  a  vertical  point  force,  Fz,  or  an  additional  contributor  to 
the  vertical  couple,  Mzz.  Recently,  Taylor  (1989)  used  the  modified  Day,  et  al. 
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model  to  argue  that  spall  is  a  significant  contributor  to  regional  Pg  and  that 
depth  of  burial  (and  hence  relative  size  of  spall)  influences  the  proposed  Mur¬ 
phy  and  Bennett  (1983)  high/low  frequency  spectral  ratio  discriminant. 

In  this  paper,  we  present  a  more  general  derivation  of  the  vertical  point 
force  at  the  surface  as  an  approximation  for  the  buried  tension  crack  represen¬ 
tation.  We  show  that  the  approximation  is  valid  in  a  more  general  context  than 
originally  shown  by  Day,  et  al.  (1983)  and  examine  the  limitations  of  the 
approximation.  Green’s  functions  for  point  source  loads  and  buried  tension 
cracks  are  examined  in  the  time  and  frequency  domain  for  a  simple  crustal 
structure.  We  see  that  the  vertical  point  force  approximation  may  suffice  for 
many  applications.  However,  the  radiation  pattern  and  frequency  dependence 
of  the  tension  crack  differs  from  the  point  force  at  high  frequencies  and  the  ten¬ 
sion  crack  representation  may  be  a  more  accurate  representation. 

Following  the  discussion  regarding  the  seismic  source  representation  of 
spall,  we  discuss  possible  smoothing  operators  to  be  used  to  represent  the 
finite  extent  of  spall.  A  take  off  angle  dependent  smoothing  operator  is  chosen 
that  has  reasonable  high  frequency  behavior,  zero  displacement  at  the  tension 
crack  tip,  and  has  a  concave  upward  spall  surface. 

Finally,  we  examine  spall  modeling  results  of  Barker  and  Day  (1990)  based 
on  axisymmetric  spall  models  that  fit  teleseismic  P  waves  radiated  by  2-D 
axisymmetric  nonlinear  finite  difference  explosion  simulations.  Regional 
seismograms  are  computed  based  on  their  spall  model  and  the  relative  excita¬ 
tion  of  Pg  and  Lg  are  compared  from  the  explosion  and  spall  sources. 
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III.  AN  EQUIVALENT  FORCE  REPRESENTATION  FOR  SPALL 

Day,  et  at.  (1983)  suggested  that  spall  may  be  represented  as  a  horizon¬ 
tally  oriented  tension  crack  that  opens  and  closes  in  the  vertical  direction  (see 
Figure  1).  Given  this  simplification,  we  can  write  that  the  far  field  displacement 
due  to  the  motion  on  a  surface  is  given  by  (Aki  and  Richards,  1980) 

UiM  =  jz  n}[ ti)  5u*(cmi)  Cjkpq(T\)  — ^ —  d\  (1) 


where  I  is  the  spall  surface  with  unit  normal,  ny(r i),  and  8uk(r\)  is  the  motion 
across  the  surface,  Cjkpq{r\)  is  the  elastic  modulus  tensor,  and  Gip(r\)  is  the 
Green's  tensor  for  the  medium.  For  the  opening  of  a  horizontal  tension  crack 
n(r i)  and  5u(r\)  are  parallel  in  the  vertical  direction,  j=z,  and  we  can  write, 


8Gfo  (ti) 

Tizz{r\)  =  C™.  (r|)  — — ,  which  by  reciprocity  is  the  vertical  stress  com- 

ponent  at  the  source  location,  rj,  due  to  a  point  force  in  the  i’th  direction  at  the 
receiver  location,  r.  Equation  (1)  becomes, 


Ui(oi,r)  =  |E  8 uz{(o,x,y)  Tizz{x,y)  dxdy. 


(2) 


We  expand  Tizz  in  a  Taylor  expansion  about  the  free  surface,  z=0, 


izz  |  z= 0 


+  z—  |  z=0  + 


dZ 


(3) 


By  virtue  of  the  free  surface,  the  first  term  is  zero,  and  Equation  (2)  becomes 

ui( co.r)  =  hs  5uz(co,x ,y)  Tizz(x,y)  \  z=0  dxdy  (4) 
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Figure  1.  The  axisymmef'C  spall  above  the  explosion  is  modeled  as  a 
circular  horizontal  tension  crack  that  opens  and  closes  in  the 
vertical  direction.  The  spall  is  parameterized  by  the  source  depth 
hs,  the  crack  radius  a,  and  the  displacement  time  function  u(t,x). 
The  radiation  from  a  horizontal  tension  crack  is  equivalent  at  low 
frequencies  to  the  radiation  from  a  vertical  point  force,  Fz,  at  the 
surface  proportional  to  the  second  time  derivative  of  5u.  The 
moment  tensor  source  representation  for  the  tension  crack  is 
proportional  to  5u.  A  is  the  area  of  the  crack,  p  is  the  density, 
and  X  and  p  are  the  Lame  parametes  for  the  medium. 


4 


where  hs  is  the  spall  depth.  Since  7}y*  is  the  stress  field  due  to  G,y,  we  may 
write  the  z’th  component  of  the  equation  of  motion, 


^  TIZ2  d  Tjyz 

d  Z  +  dy 


dT, 


1X2 


dx 


=  -co2p  Giz. 


(5) 


If  we  apply  Equation  (5)  to  the  free  surface,  z=0,  then  the  partial  derivatives 
with  respect  to  x  and  y  are  zero  since  Tiyz  =  Tixz  =  0  for  z=0,  and  we  have  that, 

— 5r- 1  m,  =  -to2pG/z(z=0).  (6) 


Substitute  Equation  (6)  into  Equation  (4), 

U/((o,r)  =  -phs032  jz  Suz(o),x,y)  Giz(x,y,z=Q)  dxdy  (7) 

Uj(w,r)  =  -p/is/4o)25u2G/2  = -msco28uzG^  = /o)/sG;2.  (8) 

Thus,  -phsco28u2(a),x,y)  is  interpreted  as  the  vertical  force  per  unit  area  at  the 
free  surface,  z=0.  Note  that  the  equivalent  force  is  applied  to  the  free  surface 
above  the  spall  surface,  and  that  the  force  has  the  time  history  of  the  spall 
acceleration.  In  contrast,  the  tension  crack  source  has  the  time  history  of  the 
spall  displacement.  The  spall  momentum,  ls,  is  given  by  the  spall  mass,  ms  = 
p hsA,  times  the  mean  velocity,  /'co5uz  and  therefore  the  equivalent  spall  force, 
Fs  =  / co/s .  In  this  manner  we  can  relate  estimates  of  spall  mass,  and  momen¬ 
tum  to  an  equivalent  point  force.  Alternatively,  we  could  use  Equation  (2)  to 
describe  the  source  as  a  buried  tension  crack, 

Uj(wj)  =  buzATizz  =  §uzA[kGix  x  +  XGjyy  +  {X+2\i)Gjz<z]  (9) 


5 


There  are  several  approximations  inherent  in  the  derivation  of  Equation  (7). 
Since  a  Taylor  expansion  was  used  from  the  free  surface,  it  is  likely  that  the 
representation  will  break  down  for  seismic  wavelengths  on  the  order  of  the  spall 
depth.  Related  to  this  limitation  is  the  fact  that  tne  radiation  pattern  for  a  point 
force  at  the  free  surface  is  not  exactly  the  same  as  that  for  a  tension  crack  at 
depth. 

We  address  these  two  limitations  by  comparing  the  regional  Green’s  func¬ 
tions  for  a  surface  point  force  and  the  buried  horizontal  tension  crack.  Equa¬ 
tions  (8)  and  (9)  imply  that  phsa£Gj2{z=0)  =  Tizz{z=hs).  Figure  2  shows  the 
Green’s  functions  for  buried  tension  cracks  compared  to  the  Green’s  function 
for  a  vertical  point  force  at  the  free  surface.  The  Green’s  functions  were  com¬ 
puted  for  the  vertical  displacement  at  a  distance  of  300  km  using  the  Eastern 
Kazakhstan  crustal  model  from  McLaughlin,  ef  al.  (1988)  (see  Table  1).  A 
wavenumber  integration  code  was  used  to  compute  the  full  response  of  the 
medium  from  0  to  5  Hz.  The  tension  crack  (Tzzz(f))  and  vertical  point  force 
(p  hs  Gzz  (f)  )  Green’s  functions  are  all  plotted  to  the  same  scale.  We  see  that 
the  relative  amplitudes  of  the  Pg,  Lg,  and  Rg  are  similar  when  we  compare  the 
point  force  and  tension  crack  representations  for  depths  less  than  500  m.  How¬ 
ever,  the  individual  wavepackets  in  the  Pg  can  be  seen  to  vary  with  depth  for 
the  tension  crack  sources. 
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VERTICAL  FORCE  VERSUS  ENSION  CRACKS  (0  TO  5  HZ 


TABLE  1. 


Eastern  Kazakhstan  Structure  Modified  from  Stevens  (1986) 


h(km) 

a(km/s) 

(3(km/s) 

p(gm/cc) 

Q a 

depth(km) 

2.000 

5.020 

2.790 

2.700 

100.0 

0.000 

1.000 

5.400 

3.000 

2.700 

150.0 

2.000 

2.488 

5.900 

3.300 

2.700 

200.0 

3.000 

10.976 

6.100 

3.400 

2.700 

600.0 

5.488 

5.488 

6.308 

3.541 

2.702 

525.9 

16.464 

5.488 

6.597 

3.703 

2.807 

500.0 

21.952 

5.488 

6.736 

3.781 

2.858 

450.0 

27.440 

5.564 

6.782 

3.807 

2.875 

400.0 

32.928 

6.504 

6.795 

3.814 

2.879 

350.0 

38.492 

8.006 

8.147 

4.573 

3.372 

179.5 

44.996 

9.359 

8.138 

4.568 

3.369 

167.5 

53.002 

10.940 

8.106 

4.550 

3.358 

159.4 

62.361 

12.780 

8.065 

4.527 

3.343 

153.8 

73.301 

14.950 

8.047 

4.517 

3.336 

150.3 

86.081 

17.470 

8.070 

4.530 

3.345 

148.7 

101.031 

20.420 

8.117 

4.556 

3.361 

148.1 

118.501 

23.880 

8.154 

4.577 

3.375 

147.9 

138.921 

27.910 

8.161 

4.581 

3.378 

147.4 

162.801 

32.630 

8.145 

4.572 

3.372 

146.8 

190.711 

38.140 

8.120 

4.558 

3.363 

146.2 

223.341 

oo 

8.101 

4.547 

3.356 

145.8 

261.481 
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Figure  3  shows  an  enlargement  of  the  Pg  waveforms  in  Figure  2.  A  visual 
comparison  shows  that  correct  phasing  occurs  between  the  point  force  and  the 
tension  crack  sources  for  low  frequencies  while  differences  are  most  pro¬ 
nounced  at  high  frequencies  and  greater  depths  due  to  interference  effects 
between  pP  and  P  components  of  the  Pg  wavetrain. 

To  illustrate  this  further,  Figure  4  compares  the  point  source  and  buried 
tension  crack  sources  in  the  frequency  domain  for  windows  taken  around  the 
Pg,  Lg,  and  Rg  phases.  We  see  that  the  shallowest  tension  crack  (100  m)  is 
well  represented  by  a  point  force  up  to  5  Hz  for  the  Pg,  Lg,  and  Rg.  Modulation 
of  the  700  m  tension  crack  source  is  clearly  present  in  the  Lg  ratio  due  to  sS+S 
interference.  It  appears  that  the  vertical  point  force  representation  is  a  good 
approximation  to  the  buried  tension  crack  at  least  up  to  spall  depths  of  1/2  the 
wavelength.  For  spall  depths  of  500  meters  this  corresponds  to  about  5.0  Hz 
for  P  waves  in  this  model  and  about  2.8  Hz  for  S  waves.  For  this  velocity  struc¬ 
ture  the  differences  at  several  slownesses  (Pg,  Lg,  and  Rg)  between  the  ten¬ 
sion  crack  and  vertical  point  force  representation  are  probably  not  significant 
below  5  Hz  for  spall  depths  shallower  than  300  m. 
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VERTICAL  FORCE  VERSUS  TENSION  CRACKS  (0  ^0  5  rZ} 


SECONDS 


An  enlargement  of  the  Pg  waveforms  of  Figure  2.  Note  that  the 
similarity  of  the  waveforms  degrades  for  depth  >  300  m. 


LOG10(T27/(Gzphw  2))  LOG10(T22/(Gzph«  2))  LOC,0(T22/(Gzphu  “)) 


RG:  TENSION  CRACK  VS  VERTICAL  POINT  FORCE 


LG:  TENSION  CRACK  VS  VERTICAL  POINT  FORCE 


PG:  TENSION  CRACK  VS  VERTICAL  POINT  FORCE 


Figure  4.  Spectral  ratios  of  the  tension  crack  seismograms  over  the  vertical 
point  force  seismogram  of  Figure  2.  Windows  were  taken  around 
the  Pg,  Lg  and  Rg  phases,  Tension  crack  seismograms  for 
source  depths  of  100,  300,  500  and  700  m  are  shown.  In 
general,  the  approximation  of  the  tension  crack  by  the  vertical 
point  force  degrades  with  depth  of  the  tension  crack. 
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IV.  THE  FINITE  EXTENT  OF  SPALL 


At  high  frequencies,  the  distributed  nature  of  spall  is  not  adequately 
represented  by  a  single  point  force  on  the  surface  or  even  by  the  moment  ten¬ 
sor  representation  of  a  tension  crack.  We  wish  to  derive  some  simple  expres¬ 
sions  for  the  high  frequency  asymptotic  behavior  of  the  spall  signal  in  the  far 
field.  Suppose  that  the  spall  displacement  is  specified  over  a  disk  of  radius,  a, 

then  for  a  wavepacket  with  horizontal  wavenumber,  k  =  to  s'n  ^  ,  where  <t>  is 

c 

the  incidence  angle,  the  far  field  displacement  at  range,  r,  will  be  approximated 
by 

Uj{a),r,k)  ~  Uj((o,r ,k)  ^  S(w,x,y)e~lkxdxdy  (10) 

3 

=  Ui(a>,r,k)  JQ  JQ  S (w,s ,Q)e~ikscos{P) sdsdQ 

Uj{(a,r,k)S((o,s,Q)  describes  the  contribution  to  the  seismogram  with  horizontal 
wavenumber  k  from  a  location  on  the  spall  disk.  If  we  assume  that 
Uj{u>,r,k)S{(Ji,s,Q)  is  independent  of  0  then  we  have  a  zero'th  order  Hankel 
transform  of  the  form, 

Uj{(a,r,k)  ~  Uj{o3,r,k)  2njQ  S{(a,s)J0{ks)sds.  (12) 

If  we  assume  that  U/(a),r,k)S(o),s,0)  is  uniform  over  the  disk  then  we  have 

u,(ou,/0  =  Ui{(a,r,k)~^-J-\(ka)  =  Ui((o,r,k)E(k),  (13) 

(ka) 

where  J^(.)  is  the  I’st  order  Bessel  function.  For  large  ka,  E(k)  is  approximated 
by  ^/2/7^a3/2(ka)■3/2cos(ka-3•J).  Therefore,  this  smoothing  function  imposes  an 
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or3/2  asymptotic  roll-off  at  high  frequencies. 


A  convenient  form  for  the  spall  distribution  is  S(co,s)  =  (a2-s2)/a2,  since  it 
has  zero  displacement  at  the  crack  tip  and  the  spall  separation  is  concave 
upward.  Under  this  assumption  the  smoothing  function  has  spectral  form, 

2 

m  =  (i4) 

{kar 

which  is  proportional  to  the  2’nd  order  Bessel  function  and  for  large  ka  rolls-off 
like  cd-5/2.  This  simple  model  is  preferred  becasue  it  has  the  physically  reason¬ 
able  property  of  zero  displacement  at  the  crack  tip,  and  the  spall  surface  is  con¬ 
cave  upward  as  suggested  by  Eisler  and  Chilton  (1964).  Because  this  continu¬ 
ous  spall  distribution  function  is  zero  at  the  crack  tip  the  smoothing  operator 
falls  off  faster  than  or2  at  high  frequencies.  Given  the  sparsity  of  information 
regarding  the  radial  dependence  of  spall,  this  simple  model  will  suffice  to  com¬ 
pare  the  relative  importance  of  the  seismic  spall  source  with  that  of  the  explo¬ 
sion  source.  It  should  be  noted  that  any  simple  continuous  distribution  function 
with  zero  at  the  crack  tip  will  fall-off  faster  than  co-2  at  high  frequencies  and 
therefore  the  far-field  signal  at  high  frequencies  should  be  dominated  by  the 
explosion  source. 
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V.  SPALL  SIGNAL  SYNTHETICS 


Barker  and  Day  (1990)  have  used  teleseismic  P  waves  predicted  by  2-D 
axisymmetric  nonlinear  finite  difference  simulations  of  explosions  to  derive 
parametric  models  for  the  equivalent  spall  source.  They  found  that  the 
differences  between  the  nonlinear  simulations  and  linear  point  explosion 
sources  could  be  modeled  by  a  simple  parametric  spall  model.  Their  model  is 
a  circular  horizontal  tension  crack  of  radius,  a,  depth,  hs,  and  has  a  uniform  ini¬ 
tial  velocity  distribution  between  v1  and  v2.  The  reader  is  referred  to  Barker 
and  Day  for  a  detailed  motivation  of  the  model.  The  function  5£7(f)  appropriate 
for  use  in  Equations  (8)  or  (9)  is  given  by 

2/4 

8 U(t)  =  t(v2+  v^)/2-gtz/2,  0  ^  t  <, — -  =  f,  (15) 
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817(f)  = 


t(v2  -  gt/2)z 
2{v2  -  v,) 


2/i  2  v2 

=  t,  <  t  £  — -  =  t2 
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The  2-D  nonlinear  axisymmetric  finite  difference  simulations  are  described 
in  McLaughlin,  et  at.  (1988)  and  are  intended  to  simulate  125  Kt  explosions  at 
200,  300,  680,  and  980  meter  depths  in  a  Shagan  River  Test  Site  structure. 
We  concentrate  on  the  680  and  the  980  meter  depth  of  burials  which 
correspond  to  roughly  normal  scaled  depth  of  burial  (122  m /Ktv3)  and  60% 
overburied  (196  m/Kf1/3)  respectively.  The  spall  parameters  determined  by 
Barker  and  Day  are  listed  below  in  Table  2.  Figure  5  shows  the  equivalent 
spall  sources  for  the  680  m  depth  simulation  as  determined  by  Barker  and  Day 
(1990)  for  teleseismic  P,  and  regional  Pn,  Pg,  and  Lg.  The  smoothing  operator 
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Pn  Spall  Source  Function  (NT) 


Pg  Spall  Source  Function  (NT) 


Figure  5.  The  equivalent  spall  source  (for  source  depth  of  680  m  from 
Table  2)  for  four  different  apparent  phase  velocities 
corresponding  to  teleseismic  P,  Pn,  Pg,  and  Lg.  The  smoothing 
operator  E(t)  simulates  the  distributed  nature  of  the  source. 
Slower  phases  are  more  sensitive  to  the  size  of  the  spall  disk. 
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of  Equation  (14)  was  used  to  mode!  the  finite  extent  of  the  spall  surface  while 
the  spall  source  time  function  was  specified  by  Equation  (15).  Apparent  veloci¬ 
ties  of  14,  8.2,  6.5,  and  3.7  km/s  were  chosen  to  represent  the  teleseismic  P, 
Pn,  Pg,  and  Lg  respectively.  Note  that  the  spall  is  characterized  by  a  high  fre¬ 
quency  opening  pulse  followed  by  a  low  frequency  unloading  and  a  small  spall 
closure  pulse.  Because  this  spall  model  opens  simultaneously  over  the  entire 
spall  disk  but  distributes  the  closure,  the  spall  closure  pulse  is  small  compared 
to  the  opening  pulse. 

The  Lg  and  Pg  spall  sources  are  smoother  and  lower  frequency  than  the 
Pn  and  teleseismic  P  spall  sources.  The  lower  the  apparent  velocity  of  the 
seismic  phase,  the  more  important  the  extended  size  of  the  source.  Spectra 
are  shown  in  Figure  6.  We  see  that  the  model  predicts  a  sharply  peaked 
source  spectrum  for  the  regional  phases  and  in  particular  the  Lg  falls-off  very 
rapidly  for  frequencies  above  2  Hz. 


TABLE  2. 

Barker  &  Day  (1990)  Spall  Models  for  Shaaan  River  Synthetics 


source  depth  (m) 

200 

300 

680 

980 

v-i  (m/s) 

1.1 

1.1 

1.1 

1.1 

v2  (m/s) 

14. 

12. 

3. 

1.3 

a  (m) 

600 

625 

650 

650 

hs  (m) 

100 

150 

200 

200 

/,  (1 012  Nt-s) 

2.31 

3.26 

1.47 

0.86 
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Figure  6.  Spectra  of  the  time  series  of  Figure  5.  Note  that  the  Lg 
equivalent  spall  source  falls-off  more  rapidly  than  the  Pg  or  Pn 
spall  source. 
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Figure  7  shows  the  tension  crack  Green’s  function  for  a  distance  of  300  km 
convolved  with  the  appropriate  Pg  and  Lg  spall  source  functions  compared  with 
the  explosion  Green’s  function  and  the  appropriate  explosion  RDP.  Note  that 
the  Lg  explosion  signal  is  insignificant  compared  to  the  Lg  spall  signal.  Based 
on  Nuttli’s  (1986)  formulas,  mb{Lg)  =  6.2  for  the  synthetic  spall+explosion  Lg 
signal  using  Q  =  350,  and  f  =  1  Hz.  In  contrast,  the  exp'osion  synthetic  mb{Lg) 
is  4.9.  Given  Nuttli’s  NTS  mb{Lg): yield  formulas  the  total  Lg  excitation  is  about 
twice  as  large  as  would  be  predicted  for  a  125  Kt  NTS  explosion. 

It  should  be  emphasized  that  the  explosion  source  is  located  in  a  high- 
velocity  (a  >  4500  m/s)  layer.  The  phase  velocities  of  Lg  are  less  than  the  P- 
wave  velocity  of  the  source  medium  and  the  explosion  is  an  inefficient  source 
for  Lg  in  such  a  model.  The  explosion  is  a  more  efficeint  source  of  Lg  for  a 
low-velocity  medium  such  as  NTS  (see  Barker,  et  at.  1990). 

In  comparison  to  the  Lg,  the  raw  displacements  of  the  Pg  spall  and  explo¬ 
sion  signals  are  of  the  same  order  of  magnitude.  Spectra  of  the  Pg  and  Lg  sig¬ 
nal  windows  are  compared  in  Figure  8.  The  spall  contribution  dominates  the  Lg 
signal  up  to  3  Hz.  At  1  Hz  the  spall  contribution  to  Lg  is  20  times  larger  than 
the  explosion  contribution.  At  1  Hz,  the  Pg  spall  contribution  is  about  twice  as 
strong  as  the  explosion  part.  The  explosion  -'ource  dominates  the  Pg  signal 
above  1 .5  Hz. 

If  we  compare  the  Barker  and  Day  spall  model  with  other  published  spall 
sources,  we  find  that  it  is  intermediate  in  strength.  Viecelli  (1973)  found  for 
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Figure  8.  Spectra  of  the  time  series  in  Figure  7.  Note  that  the  Lg  spall 
signal  is  20  times  larger  than  the  Lg  explosion  signal  at  1  Hz  and 
the  Pg  spall  signal  is  about  two  times  larger  than  the  explosion 
Pg  signal  at  1  Hz. 
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explosions  at  normal  scaled  depth  of  burial  that  spall  momentum  scaled  linearly 
with  yield  as  4.6*  109  Nt-s/Kt.  The  Barker  and  Day  spall  model  for  the  680 
meter  depth  of  burial  corresponds  to  about  2.5  times  the  spall  predicted  by 
Viecelli’s  scaling  relationship.  This  is  about  the  same  factor  found  between 
Viecelli’s  scaling  relationship  and  the  2-D  nonlinear  finite  difference 
PILEDRIVER  simulations  described  in  Day,  et  a!.  (1983).  Sobers  (1978)  scal¬ 
ing  relationships  predict  spall  momenta  several  times  larger  than  the  Barker  and 
Day  model.  It  appears  that  the  2-D  Shagan  River  nonlinear  simulations  predict 
spall  momenta  roughly  intermediate  between  the  Viecelli  (1971)  and  Sobel 
(1978)  scaling  relationships. 

Patton  (1988)  estimated  spall  momentum  for  the  Pahute  Mesa,  NTS, 
explosion  HARZER,  to  lie  between  1  and  2  *  1012  Nt-s  with  an  explosion 
moment  between  7  and  5.7  *  1015  Nt-m.  The  simulation  shown  in  Figures  7 
and  8  had  an  explosion  moment  of  about  1.4  *  1016  Nt-m.  The  spall  momen¬ 
tum  for  the  simulation  is  therefore  a  factor  of  about  two  lower  than  Patton’s  esti¬ 
mate  *cr  HARZER  based  on  simple  linear  moment  scaling. 

Figures  9  and  10  show  the  Lg  and  Pg  waveforms  and  their  spectra  for  the 
980  meter  depth  of  burial  simulation  using  the  spall  parameters  listed  in  Table 
2.  Note  that  in  comparison  to  the  680  meter  depth  of  burial  calculation,  the  980 
meter  deep  explosion  contribution  to  Pg  is  larger  than  the  spall  contribution. 
The  spall  produces  only  5  to  7  times  more  Lg  at  1  Hz  than  the  explosion,  and 
the  bandwidth  is  narrower  for  which  spall  is  a  greater  source  of  Lg  than  the 
explosion.  The  total  Lg  amplitude  at  1  Hz  for  the  980  meter  depth  of  burial  is 
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one-fourth  the  Lg  amplitude  for  the  680  meter  simulation. 
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Lg  WAVEFORMS  (H=980  M,  A  =  300  KM) 


SPALL+EXPL 


SECONDS 


Pg  WAVEFORMS  (H=980  M,  A=300  KM) 


SPALL+EXPL 


SECONDS 


Synthetic  Lg  and  Pg  waveforms  at  a  distance  of  300  km  for  the 
980  meter  depth  of  burial  1 25  KT  simulation.  The  Lg  spall  signal 
dominates  the  Lg  explosion  signal  while  explosion  Pg  signal  is 
slightly  larger  than  the  spall  Pg  signal. 
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FREQUENCY  (HZ) 

Figure  10.  Spectra  of  the  time  series  in  Figure  9.  Note  that  the  Lg  spall 
signal  is  about  5  to  7  times  larger  than  the  Lg  explosion  signal, 
while  the  explosion  Pg  signal  is  larger  than  the  spall  signal  at 
nearly  all  frequencies. 
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DISCUSSION  AND  CONCLUSIONS 


We  have  presented  a  more  general  derivation  that  a  vertical  point  force  at 
the  surface  may  approximate  the  source  of  seismic  radiation  for  a  buried  ten¬ 
sion  crack.  The  approximation  was  tested  for  a  simple  velocity  structure,  and  it 
was  found  valid  for  frequencies  up  to  5  Hz,  and  spall  depths  shallower  than  500 
meters.  One  caveat  to  this  conclusion  is  that  the  Green’s  functions  were  com¬ 
puted  for  a  fairly  high  velocity  structure  with  a  thick  surface  layer.  The  frequen¬ 
cies  for  which  the  approximation  is  valid  may  be  lower  for  lower  velocity  source 
structures  or  more  complicated  velocity  structures. 

A  smoothing  operator  that  simulates  the  distributed  nature  of  the  spall  was 
presented.  The  smoothing  operator  (filter)  has  the  desirable  characteristics  that 
it  has  zero  motion  at  the  crack  tip,  and  the  spallation  surface  is  concave 
upward,  therefore  the  spectrum  falls-off  at  high  frequencies  faster  than  or2.  A 
specific  model  is  presented  that  falls-off  at  high  frequencies  proportional  to 
o r5/2  Because  the  spall  signal  is  proportional  to  to2  for  low  frequencies,  the 
spall  signal  is  expected  to  be  a  narrow  band  signal  compared  to  the  main 
explosion  signal. 

The  spall  model  and  parameters  of  Barker  and  Day  (1990)  were  used  to 
synthesize  regional  waveforms  for  Pg  and  Lg  and  the  spall  contribution  to  the 
seismogram  was  compared  to  the  point  explosion  contribution.  The  Lg  and  Pg 
signals  were  found  to  be  significantly  affected  by  the  spall  source  for  the  2-D 
axisymmetric  nonlinear  finite  difference  simulations. 
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We  reach  a  conclusion  similar  to  Taylor  (1989)  that  spall  is  probably  a 
significant  contributor  to  the  regional  Pg  signal.  This  is  not  unexpected  since 
spall  partially  cancels  the  pP  and  replaces  it  with  a  lower  frequency  and 
delayed  spall  signal.  In  addition  we  conclude  that  spall  is  almost  certainly  a 
significant  contributor  to  the  Lg  signal  for  high-velocity  crustal  structures.  It  is 
important  to  note  that  for  explosions  in  a  low-velocity  (a  <  4500  m/s)  near¬ 
surface  layer,  the  explosion  excitation  of  Lg  is  more  efficient.  The  relative  exci¬ 
tation,  explosion  versus  spall,  of  Lg  is  not  as  divergent  in  a  low-velocity  source 
structure  as  in  a  high-velocity  source  structure.  These  conclusions  require 
adjustment  for  NTS  structures  (see  Barker,  et  al.  1990). 

It  should  be  noted  that  the  spall  model  used  in  this  paper  contains  approxi¬ 
mations  and  lacks  characteristics  that  should  become  important  at  higher  fre¬ 
quencies  (>5  Hz)  than  were  examined  in  this  paper.  First  of  all,  the  spallation 
does  not  occur  across  a  simple  surface.  The  spall  is  distributed  over  a  volume 
and  therefore  the  spall  parameters  ( v-[,v2,hs,a )  represent  integrals  over  a  non¬ 
linear  volume.  Caution  should  be  used  in  making  too  detailed  physical  interpre¬ 
tations  of  these  model  parameters.  Second,  the  entire  spall  surface  does  not 
open  simultaneously.  This  will  further  smooth  the  initial  spallation  signal.  Third, 
the  spall  is  not  axisymmetric  and  it  may  radiate  SH  waves  as  well  as  have  non 
axisymmetric  radiation  patterns  for  P  and  SV  waves. 

It  is  clear  from  our  simulations  that  the  spall  source  will  depend  on  scaled 
depth  of  burial.  Overburied  explosions  will  tend  to  spall  less  mass  for  a  shorter 
time  and  hence  have  a  smaller  amplitude  and  higher  frequency  spall  signature. 
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We  should  expect  that  spall  scaling  will  also  depend  on  material  properties  and 
structure  above  the  emplacement  as  well  as  the  properties  of  the  rock  near  the 
working  point.  Since  the  effects  of  gravitational  forces  on  containment  do  not 
scale  linearly  with  yield,  we  also  expect  that  spallation  may  be  more  important 
for  smaller  explosions.  For  example,  Stump  (1985)  found  for  chemical  explo¬ 
sions  in  alluvium,  that  the  spall  mass  was  2  to  3  times  larger  than  that  predicted 
by  Sobel’s  scaling  relationship.  This  may  have  been  due  to  either  weaker 
material  properties  or  lack  of  linear  scaling  of  gravitational  effects  with  yield. 

The  mb( Lg)  predicted  for  the  normal  buried  simulation  (Figure  7)  is  roughly 
0.6  to  0.7  magnitude  units  larger  than  the  mb(Lg)  predicted  for  the  60%  over 
buried  simulation  (Figure  8).  If  spall  is  the  dominant  source  for  Lg  excitation 
then  we  should  expect  mb( Lg)  to  depend  upon  scaled  depth  of  burial.  However, 
given  the  robust  nature  of  mb(Lg):yield,  it  is  likely  that  other  mechanisms  are 
responsible  for  at  least  some  of  the  Lg  excitation. 

The  presence  of  transverse  Lg  from  explosions  implies  that  non  axisym- 
metric  nonlinear  mechanisms  or  scattering  are  also  responsible  for  Lg  excita¬ 
tion.  Although  non  axisymmetric  spall  may  excite  transverse  Lg  (Lg-SH),  it  has 
been  argued  by  Gupta  and  Blandford  (1983)  that  the  Lg-SH  signal  grows  rela¬ 
tive  to  the  vertical  and  radial  components  of  Lg  (Lg-SV)  out  to  a  coupie  of  hun¬ 
dred  km.  This  may  imply  that  scattering  is  operating  along  the  path  but  does 
not  in  itself  explain  the  initial  excitation  of  Lg-SV.  Conversion  of  compressional 
energy  may  be  converted  at  interfaces  with  significant  velocity  contrasts  besides 
the  free  surface.  Or,  scattering  by  near  source  lateral  heterogeneity  may 
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contribute  to  the  excitation  of  Lg  and  Lg-coda. 

A  number  of  overburied  explosions  have  excited  Lg  as  efficiently  as  nor¬ 
mally  buried  explosions.  These  explosions  have  been  located  in  diverse  geolo¬ 
gic  environments  at  scaled  depths  of  burial  greater  than  normal  containment 
depth  (~120  m/Ktv 3)  and  include  SCOTCH  (180  m/Ktv 3  in  tuff),  REX  (250 
m/KtV3  in  tuff),  SALMON  (475  m/KtV3  in  salt),  GNOME  (250  m/Ktv 3  in  salt), 
RULISON  (750  m/Ktv 3  in  shale),  GASBUGGY  (420  m/KtV3  in  shale),  and  RIO 
BLANCO  (>420  m/Ktv 3  in  sandstone)  (yields  and  depths  from  Springer  and 
Kinnaman,  1971  and  1975).  Nuttli  (1986)  has  estimated  the  yields  of  six  of 
these  events  to  within  50%  based  on  mb(Lg),  and  Blandford,  et  al.  (1981) 
have  shown  that  the  relative  Lg  excitation  from  GNOME  and  SALMON  is  con¬ 
sistent  with  their  yield  ratio.  The  Lg  excitation  by  these  seven  events  argues 
that  there  are  other  mechanisms  that  contribute  to  Lg  since  spall  was  reduced 
for  these  explosions.  A  number  of  these  events  were  located  in  geologically 
complicated  environments  and  scattering  may  have  been  important  (SALMON, 
GNOME)  for  some  of  the  events.  Events  at  NTS  (SCOTCH  and  REX)  were 
detonated  in  low-velocity  materials  (a  <  4500  m/s)  and  the  relative  Lg  contribu¬ 
tion  of  spall  and  explosion  is  more  nearly  equal  in  an  NTS  structure.  RIO 
BLANCO  consisted  of  three  30  Kt  explosions  separated  by  1 20  m  and  therefore 
the  nonlinear  zone  may  have  deviated  from  spherical  symmetry  which  would 
have  introduced  a  CLVD  component  to  the  source  and  excited  Lg.  Furthermore, 
moderate  spall  was  documented  for  RIO  BLANCO  by  Toman  et  al.  (1973). 
Murphy  and  Archambeau  (1986)  argue  that  RULISON  was  contaminated  by 
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tectonic  release  in  the  short-period  frequency  band  and  the  Lg  excitation  may 
reflect  this  deviatoric  source.  Murphy  and  Archambeau  also  found  it  necessary 
to  incorporate  spall  into  the  modeling  of  teleseismic  P  waves  for  GASBUGGY 
and  RULISON  even  though  these  events  were  overburied.  In  addition,  GAS- 
BUGGY  and  RULISON  were  detonated  in  layered  material  with  P-wave  veloci¬ 
ties  less  than  or  equal  to  4500  m/s  so  the  explosion  source  could  excite  Lg 
directly.  Consequently,  there  appear  to  be  no  clear  cut  cases  where  an  over¬ 
buried  explosion,  uncontaminated  by  non-spherical  effects,  was  detonated  in 
high-velocity  layered  rock  for  which  we  have  mb{Lg )  measurements. 

Based  on  nonlinear  finite  difference  simulations  for  explosions  in  brittle 
high-velocity  rock  we  conclude  that  spall  should  be  a  significant  source  of  Lg. 
Other  mechanisms  such  as  near-source  scattering  may  also  be  viable  mechan¬ 
isms  for  Lg  excitation.  A  theoretical  basis  for  yield  estimation  or  discrimination 
based  on  Lg  will  require  a  better  understanding  of  these  mechanisms. 
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